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are well known, the relative ease with which the Casimir force density in even complex 
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herein. Results may be seen as an important building block in the settling of issues 
of fundamental interest, such as the long-standing dispute over the thermal behaviour 
of the Casimir force or the question of what is the correct stress tensor to apply, a 
discussion re-quickened by the newly suggested alternative theory due to Raabe and 
Welsch. 
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1. Introduction 

Over the last decade or so enormous progress has been made in experimental tests of 
the Casimir effect 0. This macroscopic manifestation of quantum electrodynamics, 
once something of a curiosity subject mainly to the scrutiny of a few theorists, has 
been measured with high precision and is now spoken of as possibly exploitable in 
nanoelectromechanical applications t 2j. 

Although only a single experiment has so far employed parallel plates and with 
moderate accuracy [3], experiments employing other geometries (typically a sphere 
and a plate) have normally had to resort to planar geometries for theoretical support, 
accompanied by the neccessary corrections to account for curved surfaces (see e.g. jl]). 

The individual pieces of theory assembled in this paper are not in themselves new; 
the paper draws heavily on several references, many of which more than a decade old. 
The theory of Green's functions in a dielectric multilayer was treated by Tomas eleven 
years ago [H] , building in turn on previous work by Mills and Maradudin [Hj two decades 
earlier. Companioned by the now classical theory by Lifshitz and co-workers j7j and 
standard optical theory of reflection it provides all the necessary tools. Desipte this 
fact however, the ease with which the Casimir forces in plane parallel systems may 
be expessed appears not to be commonly recognised, although it has been implicitly 
employed by Tomas ([B] and later papers). Furthermore, no publication exists to the 
author's knowledge, explaining explicitly the procedure derived and demonstated herein. 

This paper provides background theory to aid the settlement of at least two ongoing 
disputes in the Casimir branch. Firstly, the as yet unsettled disagreement over the 
temperature effect of the Casimir force (for a recent discussion see [9J and references 
therein); numerical and theoretical treatment of the expressions obtained e.g. in chapter 
HJusing dispersion data for real materials provide predictions to settle experimentally the 
existence or non-existence of the large thermal variations of the force upheld by many. 
Secondly, doubts have been raised recently about the applicability of the Minkowski 
stress tensor and an alternative, Lorentz force-based tensor was suggested [TP] , in 
turn disputed by Pitaevskii [TT]. A procedure similar to that presented here seems 
to have been employed by Tomas in his calculations of the effects of the Lorentz-type 
tensors ^21 EH] ■ The discussion of the applicability or indeed correctness of their theory, 
however, is not within the scope of this paper. 

We have structured the paper as follows. In chapters El to 0] the background theory 
of Green's function calculation of the Casimir attraction is derived briefly, arriving at 
equation l)15J). the main result of the paper. In chapter El we demonstrate the strength 
of the procedure by using it to readily reproduce an array of previous results in various 
configurations: two half-spaces, a plate and a wall, a plate in a cavity and two plates. 

Many detailed calculations which are straightforward in principle have been 
omitted. For details, the reader may refer to [13] . 
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2. Background theory: force on an interface 

When electro- and magnetostrictive contributions are neglected, forces acting inside 
magnetodielectric media assuming no net external charge or currents are present, may 
in general be expressed through space components of the Minkowski energy momentum 
tensor Assuming isotropic, homogeneous and linear media the electromagnetic 

force density acting at position r is 

fi(r) = d k T ik = -U^d^r) - ^ H 2 d lf i(r), (1) 

where T ik is the Maxwell stress tensor, 

Tik = EiD k + HiB k — ^5 ik (E D + H B), (2) 

where indices i, k e {x, y, z) denote Cartesian vector components. We have suppressed 
the frequency dependence of the permittivity and permeability, respectively e and //, 
both defined relative to vacuum so that Di = e^eEi and Bi = fi p,Hi. 

We introduce the classical Green's dyadic T ik according to the convention of 
Schwinger and co-workers ^H] defined according to 

E(x) = — f dV F (x, x') ■ P(x') (3) 
eo J 

where x = (r,t). Due to causality, t' is only integrated over the region if < t: the 
polarisation at a time t' cannot influence the resulting electric field at time t prior to t' . 
The definition © ensures that T is a generalised susceptibility. It is well known that 
according to Maxwell's equations, V satisfies 

V x Vx T (r, r>- u) - ? (r , r ' ; u) = ^K 5(r _ r >) T, (4 ) 

c c 

where we have performed a Fourier transformation according to 

T(x,x')= r%-+» ?(r,r», (5) 

J-oo ^ 

with r = t — if. 

Invoking the fluctuation-dissipation theorem in a standard manner (e.g. [7j), yields§ 
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i(Ei(r)E k (r')) u = - coth ( 7T7 -^J 3m {T lk {r, r'; u)} (6a) 

-) -^Curl^Curl'fc^m {T^r, r'; u)} , (Qb) 

where we have used the notation Curl^ = e^dj (e^ being the Levi-Civita symbol and 
summation over identical indices is implied), Curl^. = e^d^ where d'j is differentiation 
with respect to component j of r', and p! = The brackets denote the mean value 

of the uj fourier component of the field component products with respect to fluctuations. 
For now let us assume T = for simplicity. 

§ Compared to ref. J? , T = -^D. 
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Following the Lifshitz procedure we introduce the useful quantities 

rf fc (r,r» = r ife (r,r»; (7a) 



c 



2 



T ik( r i r '\ = -z Curl ii Curl' fcm r im (r, r'\ u), (7b) 



u 2 

and replace the (fluctuation averaged) field components contained in with their 
Green's function equivalents. 

Consider now a sharp interface between two different magnetodielectric media 
such as described above. The force density acting on the interface equals the net zz- 
component of T ik when r — > r' at the surface, averaged with respect to fluctuations, that 
is the net flow of momentum through the interface when all unphysical contributions 
have been subtracted, i.e. all that flowing in equal amounts both ways across the 
interface. The resulting force expression is [7] (the original expression has been 
generalised to allow /i ^ 1) 



^\ xx yy zz I ' \ xx yy zz 



o \ xx 1 yy zz / 1 . \ xx 1 yy zz ^ , • (8) 

where a standard frequency rotation to uj = i( has been performed. The dependence 
of the components of F on r, r' and u has been suppressed and superscript h indicates 
that only the homogeneous solution of (UJ) is included. A solution of (UJ) generally takes 
the form f(z — z') + g(z + z') and the superscript ' — ' denotes that terms of T dependent 
on z + z' are discarded! . The limit r — > r' is taken so that (JHJ) is evaluated entirely on 
one side of the interface. 

Prior to solving (J3J) explicitly in some planar geometry, we introduce one further 
Fourier transformation: 

? (r,r» = J ||Le^-(^-l) ? (z, z';k x ,u), (9) 

where the subscript J_ denotes a direction perpendicular to the longitudinal or z-axis, 
i.e. in the xy-plane. 



3. Green's functions in a multilayer geometry 

We regard a multilayered geometry with a total of n + 1 layers such as depicted in 
figure n Each layer is assumed homogeneous, isotropic and of infinite transverse size. 

|| The terms of T omitted from (JSJ) are non-physical in our formalism. The particular solution of 
|J1| equals the Green's function in an infinitely large and homogeneous magnetodielectric, and is thus 
geometry independent. The z + z'-terms, argues Lifshitz, make no contribution to the net flux of 
momentum inside a single homogeneous medium, and it is straightforward to show formally that they 
make zero contribution to (JSJ) if included ^3] . These terms may be understood as a constant background 
of radiation contributing to what Milton in his book refers to as 'bulk energy' JSJ, generally of different 
value in different materials. Notably, if one followed the procedure Schwinger, Milton et.al. use to 
include the stress tensor, all terms of the homogeneous solution must be included since the Green's 
function is evaluated on both sides of an interface, i.e. in two different media. Both approaches yield 
the same end result as they should. 
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Figure 1. The multilayer geometry as discussed. We need to calculate the 
homogeneous part of the Green's function in layer (j). 

The thickness of some layer (I) is denoted a;. The Casimir attractive force per unit 
transverse area between multilayers (+) and (— ) bordering on either side of some (not 
arbitrarily chosen) layer (j) so that < j < n may be calculated using (jHJ by evaluating 
the Green's function near any of the boundaries of layer (j). In a real setting, layers 
(0) and (n) will both typically be vacuum or air, and it should be obvious that the two 
multilayers exert forces on each other in a resiprocal manner (there is a subtlety when 
material does not equal material n as discussed below). 

To mas shows jH| how the homogeneous solution of (JH) in layer (j) is found in the 
(kj_, uj) Fourier domain (now treating z and z' as parameters and k± and uj as variables, 
not the other way around) ^f. We have repeated his procedure allowing /i(r) ^ 1, yielding 

h 2 ™ -n-a- 

q=TE W 

Here z, z' G (j) and we use the notation AB = A® B. The vectors 

£>{k u u>; z) ee e+(k x )e-"^ + r+e^(k ± )e K ^\ (10) 
£<(k ± , u; z) = e~(k ± )e^ + r~e+ (k ± )e~^ (ll) 

describe waves propagating in layer (j) towards the right and left respectively and which 
are reflected off the bordering interfaces. Here 

e| Ej (fc ± ) = k± x £ = -e^ Ej (-k±) 
If For the sake of comparison, ref. |S] makes use of the quantity Go- = ^f-r^- 
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are direction vectors in a cartesian coordinate system (k±, z, kj_ x z) defined relative to 
the wave vector: k = k±k^_ + knz. Futhermore one has introduced the quantities 



k 



3 

93 ~ x ' qj' qj y 



£g — 5g,TM — 5 9j TE- 

For physical reasons, recognising that k\\( z )\ ze ^ — 1K j we must choose Shn{«j} < 0. 
The polarisation mode q runs over the polarisations TE and TM and the Fresnel 
coefficients r ■ are the relative reflected amplitudes of a g-polarised field from the entire 
stack of layers to the left (— ) or right (+) of layer (J). 

We may thus write the attraction between multilayers on either side of (J) as (it is 
sufficient to regard an interface on one side of (j) since attraction is resiprocal) 



d( f d 2 k ± /pB,h- , -pE,h- _ pB,h-\ _|_ }^(-nH,h- , -pH,h- _ -pH,h-s /, n\ 
e \ xx J ' yy,j zz,j ) ' \ xx,j ' yy,j zz,j ) \ ) 



(27T) 

where the Green's function components are taken in the limit z — > z' e (J) close to 
either of the interfaces bounding on (j). We assume in this expression that a,j may be 
varied whereas the thicknesses of other layers are treated as parameters. 

The rather complicated expression for T^ k ■ above may be vastly simplified for our 
purposes. We introduce ordinary co-ordinates according to the convention of Schwinger 
et.al. jinj by choosing x = k± so that (k±,z,k± x z) — ► (x,z,—y). Introducing the 
important quantity 

and keeping only terms dependent on z — z' we show with some lengthy but 
straightforward manipulation that the Green's function may be written very elegantly 
as 

+ h ~ r k\ 1 Kj 1 fljUJ 2 1 



IV (k ± ,uj;z,z') 



, zz — i — xx H — 21 — yy 

_€jKj OtMj £j UTMj KjC «TEj 



cosh(z — z') 



+— — (zx + xz) sinh(z — z). (14) 

It is now simple matter to calculate Yf k and Tf k and take the limit z — > z' to find 
that our final result at zero temperature becomes beautifully simple: 

j- poo poo -. 

<^°K)> = - A / dC / dk ± ■ k ± n 3 £ — . (15) 

JO JO «=TE M 

At finite temperatures, the integral over all positive imaginary frequencies in (|15|) 
becomes the sum of the residues of the coth-factors of (J6~ajl and (|6]jJ), the Matsubara 
frequencies i( m = 2nik B Tm/h, 

, T °c ' oo tm 

(^(ai)> = -— E / ^-^Er ( 16 ) 

rn=0 J{i <j=TE 93 
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Figure 2. Multiple reflections between interfaces. The scrambled area represents all 
possible paths originating and ending in layer (i). 



where the prime on the summation indicates that the m = term is given half weight. 
4. Generalised reflection coefficients 

The task remaining is to evaluate the generalised Fresnel reflection coefficients of a 
stack of magnetodielectric layers. 

For a single interface between media (i) and (j), the reflected amplitude ratio of 
a wave arriving from (i) and is partly reflected back into (i) is found from Maxwell's 
equations to be [17. 

a K i ~ lg,ii K i (17) 
r q ,% 3 - ^q,t 3 K , K LI ) 

with k as defined above and 

j fii/fij] q = TE 
\ et/ef, q = TM 

Secondly, the reflection coefficient of a system of two interfaces such as depicted in 
figure El may be calculated as the sum of coefficients pertaining to each of the infinitely 
many optical paths originating and ending in (i). We let the transmission coefficient of 
a wave transmitted from (i) to (j) be tij (omitting for now polarisation q) and recognise 
the longitudinal wave vector in (j) to be fcy = ^EjUjU) 2 / c 2 — k\ = iKj so that a wave 
travelling a distance a,j must be multiplied by exp(iA;yaj) = exp(— /%Oj); a phase shift if 
Kj is imaginary (propagating wave) or attenuation if Kj is real (evanescent wave). Thus: 



n=0 



r,-,- + r ik .e 2K i a i , . 

= / J J — ^ 18 

1 + r^r^e ZK i a i 

where we have made use of the properties t^tji — r^r^ = 1 and = —r^. Equation 
(j!8|) is valid for either polarisation respectively. 

This provides a simple procedure for calculating the reflection coefficient of a 
multilayer containing any finite number of interfaces. To calculate r~, say, as it appears 
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Figure 3. The trizone configuration of two half-spaces separated by a gap region. 



in figure we start with the leftmost interface between zones (1) and (0) and find 
r 9j io = A 9i io and invoke (fTHj) recursively to find the reflection coefficient of the two 
leftmost interfaces, then the three leftmost and so on until the closest interface, between 
(j) and (j — 1), is reached. 



5. The Casimir attraction calculated in various configurations 

We go on to demonstrate the strength of the above procedure by calculating the Casimir 
force in an array of different plane parallel configurations so as to reproduce the results 
of several references. 



5.1. Two half-spaces 

Consider first the simplest system of two half-spaces of some magnetodielectric material 
separated by a gap of width a, generally made of some other material. We denote the 
half-spaces 1 and 2 (see figure EJ) and the gap 3. 

There is now only one interface on either side of the gap and obtaining the force 
expression is almost trivial. The reflection coefficients to the right and left are 

r+ = A 9i3 2 and r~ = A ?;31 

as defined in (fTT|) . We get the force density expression (suppressing the averaging 
notation henceforth) 



fc poo j* oo ^ lv± A A „-2K 3 a 

J=°(a) = ~f d( dk ± -k ± n 3 y T \ , , as) 

V) 2 Wo Vo 3 ^ E l-A gi32 A gi3 ie- 2 ^' 1 } 

which is the classical Lifshitz result quoted in numerous references, e.g. |?1 116j. 



5.2. A plate outside a wall 

We consider a plate of finite thickness b and 
a from an infinitely thick wall of material 1 
recently been considered by Tomas ^Hj for 



material denoted 2 separated by a distance 
as depicted in figure HI Such a system has 
the sake of discussing the consequences of 
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Figure 4. The four-layered configuration of a plate and a wall. 



using an alternative Lorentz force stress tensor, and his results using the Minkowski 
tensor agrees with ours. The gap material, denoted with subscript g is allowed to be 
different from the material of the exterior, denoted with subscript e. The Green's 
function is calculated in the gap region, and r~ and r+ are the relative reflected 
amplitude of a wave originating in the gap and propagating towards the left and right 
respectively. 

We see immediately that r+ = (k 9 — 7g )ff i/ci)/(« ff + 7g, fl i«i) = — Aq,i g and combining 
(HZ!) and HE) readily yields 

_ = -A g , 2g + A g , 2e e- 2K2b 
r i 1 - A q , 2g A q , 2e e~^ b ' 

so that 

J_ = (A g , lg A g , 2g - A g , lg A g , 2e e- 2 ^ fc ) e' 2 ^ 

d qj 1 - A g , l9 A gi2s e- 2 ^ b - ( A q>lg A g , 2g - A g , l9 A g , 2e e- 2 ^ fe ) e" 2 <^ ' 

from which the Green's function and Casimir attraction per unit area between plate 

and wall follow neatly from ()14|) and (|15j). 

5.5. ^4 slab in a cavity 

We go on to study the five-layered system of a slab between two walls as defined in figure 
El A general system of five layers was first considered by Zhou and Spruch [TH] and the 
special case of a slab in a cavity was treated, apparently using the method presented 
here, by Tomas S\. Some manipulation shows that our result co-incide with the finds 
of both references. We calculate the attraction with respect to each of the gaps in turn 
and subsequently find the force acting on the slab as the difference between these. The 
procedure is identical to that above so the details of the somewhat more lengthy but 
principally uncomplicated calculations are left out. 

We denote the left and right hand gaps with superscript + and — respectively for 
reasons which will become obvious and for simplicity we assume both walls to be made 
of the same material, denoted 1, whereas the slab is made of a material indexed 2. 
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Figure 5. The system of a slab in a cavity. 



Furthermore, we use the simplifying notation A iq = A 9ji9 with q = 1, 2. With this we 
find for the left and right hand gaps 

1 UT e - 2K v a± 



where 



Uf = A lq A 2q (l - A lq A 2q e- 2 ^ a± ) - A lq (A 2q - A lq e- 2 ^ a± ) e - 2 ^\ 

r ± — 1 _ A. A. - 2 «9 a± _ A. (A„ _ A. -2K„a±\ -2 K2 b 



= 1 - A l9 A 29 e- 2 ^ - A 2q (A 2q - A lq e- 2 ^ a )e 

and resulting force on the centre slab is accordingly (a function of either a + or a~ when 
b and c are assumed constant parameters) 

Notice how, if we let the gap be either very wide (b — > oo) or perfectly reflecting 
(«2 — > oo), the terms contaning the factor exp(— 2n 2 b) vanish, and we get back Lifshitz' 
expression for two separate gaps as we should. 

A rather more instructive expression is obtained if the position of the slab is given 
not by a + and a~, but the deviation 5 of the centre of the slab from the midline of the 
cavity. Introducing the quantity h = c — b = a + + a~, we write a 1 * 1 = h/2 ± 5 and the 
force density expression becomes after some shuffling of symbols 

h f 00 . . r° . ^ A,sinh2/€„<5 



with 



q=TE 



B q = l- A 2 2q e- 2K2b + Alq(Aj q - e -^ b ) e ~^ 9 h . 



A q = 2A lq A 2q (l-e- 2 ^ b )e-^ h , 
S 2q e-^ b + A 2 lq (A 2 2q 
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Figure 6. A system of two plates of finite thickness. 



One should note in this context that this is the geometry treated by Raabe and 
Welsch [10J and later by Tomas [12] . making use of their alternative, Lorentz-type stress 
tensor. 

5.4- Two plates of finite thickness 

The system of two plates each of finite thickness has been treated by numerous authors. 
First to do so was Kupiszewska ^H] , who employed an effectively one-dimensional model 
by insisting that waves be reflected at normal incidence, as have several authors after 
her. A three dimensional geometry, allowing nonzero values of k±, appears first to have 
been considered by Jaekel and Reynaud in 1991 [20J whose result is found to be smaller 
than ours by a factor 1/2 (it should be noted that comparison is not trivial due to formal 
differences). A number of other references [TS] I2*T] 12*2"]. however, obtain results agreeing 
perfectly with that found by using the above procedure. The principally identical system 
of two semispaces each covered with a thin layer of a different substance was considered 
by Klimchitskaya and co-workers j2S], again in agreement with the below. 

The geometry is defined in figure |B] Just as before we determine the reflection 
coefficients of either slab as seen from the gap, 

+ = -A g , 2g + A g , 2e e- 2K2b2 _ -A qAg + A t?;le e- 2Klbl 

r * 1 - A 9i2s A g , 2e e- 2 ^ r<? 1 - A g , lg A 9 , le e-^ ' 1 j 

and the Casimir attractive force per unit area of the plates is found neatly from (|13j) 
and (HSJ). 

6. The neglected bulk force 

In the case that materials and n in figure ^ are different, a second force appears in 
addition to the Casimir force. It is a force acting on the entire multilayer system and 
we shall refer to it as a bulk force. Quantitatively the electromagnetic force density was 
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found to be / =T • V and the force acting on some volume V is 
F = I d 3 rf = I T -dS 

Jv JdV 

where the divergence theorem has been invoked and dS points normally out of V, 
enclosed by the surface dV. Let V be a box enclosing all interfaces of the multilayer 
system so that its ^-boundaries lie at z~ e (0) and z + e (n). Only the sides of V parallel 
to the xy-plane contribute to the bulk force, which is evaulated per unit transverse area 
as (averaging with respect to fluctuations is understood) 



bulk 



t zz (z + : 
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z/j 
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El -El 
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If now media (0) and (n) are the same, the mean squared fluctuating fields will be 
the same on either side of the multilayer, and the bulk force is zero, otherwise J^buik is 
generally nonzero. 

This force is typically neglected, and an argument in favour of doing so is surely 
that in a real system, layers are not infinitely thick. The outmost layers of a multislab 
configuration should realistically be air or vacuum, and if they are not, it simply means 
that V does not contain the entire system and the bulk force is identically cancelled by 
reflections at surfaces (not necessarily parallel with the system, or even plane) outside 
V, as it should be according to Newton's third law. 



7. Conclusion 



With the above procedure, the calculation of Casimir forces in even complex multilayered 
geometries is both quick and straightforward and shown able to reproduce the results of 
a number of previous works. Various configurations may thus be considered theoretically 
and numerically with ease to study the various dependencies e.g. on material properties 
and temperatures. The procedure may furthermore be repeated to reveal differences 
between various electromagnetic stress tensors, a subject of dispute for decades [13 121] • 
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